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Abstract 

We investigate spinodal decomposition kinetics of an initially nematic dispersion of rod-like 
viruses (fd virus). Quench experiments are performed from a flow-stabilized homogeneous nematic 
state at high shear rate into the two-phase isotropic-nematic coexistence region at zero shear rate. 
We present experimental evidence that spinodal decomposition is driven by orientational diffusion, 
in accordance with a very recent theory. 
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I. INTRODUCTION 



Systems that are quenched into a state where at least one order parameter is unstable un- 
dergo spinodal phase separation. Here, the initially homogeneous system is unstable against 
fluctuations of arbitrary small amplitude, and phase separation sets in immediately after 
a quench. In the initial stage of phase separation an interconnected "labyrinth structure" 
of regions with somewhat higher and lower values of the order parameter is observed. For 
systems containing spherical particles the relevant order parameter is the concentration. 
As Onsager showed in 1949 1], the situation is different when the particles are not spher- 
ical in shape, i.e. disk-like or elongated particles. Here the system can become unstable 
or meta-stable with respect to fluctuations in orientation. These orientational fluctuations 
drive concentrations differences, resulting in a phase with high concentration and orienta- 
tional order, the nematic phase, and a phase with low concentration and no orientational 
order, the isotropic phase. For very long and thin rods with short-ranged repulsive inter- 
actions, the binodal concentrations, i.e. the concentrations of the isotropic and nematic 
phases in equilibrium after phase separation is completed, have been determined using dif- 
ferent approximations in minimizing Onsager's functional for the free energy (see ref. Q] 
and references therein) , while for shorter rods computer simulations have been performed to 
obtain binodal concentrations jsH^J. The_spinodal concentration, where the isotropic phase 
becomes unstable has been obtained 
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Recently, a microscopic theory was developed by one of the authors, describing the initial 
stage of the kinetics of spinodal decomposition (SD) 6]. It is shown there that demixing is 
dominated by rotational diffusion and not by translational diffusion as suggested in earlier 
work [7, This is in line with results based on Ginzburg-Landau equations of motion, 
where the importance of the coupling between concentration and orientation was studied 

Q] for rod-polymer mixtures. 

Signatures of SD have been observed for suspensions of boehmite rods, by homogenizing 
a phase separated system and se quen tial polarization microscopy and Small Angle Light 
Scattering (SALS) measurements [11]. For such experiments, however, there is always an 
experimental lapse time between homogenization and the first moment of observation. In 
these experiments the initial state is not well defined. Ideally one would like to perform a 
concentration quench from low or high concentration into the two-phase region, where the 
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initial state is isotropic or nematic, respectively. In a recent paper such a kind of 'quench 



was performed by inducing polymerization of short actin chains 12j . Alternatively, external 
fields like shear flow or magnetic fields [1^ can be used to achieve well-defined quenches. 
Switching on or turning off such an external field can take the initially homogeneous system 
into either an unstable or meta-stable state. 

In this paper we induce a nematic phase with a well defined director by imposing shear 
flow to a dispersion of colloidal rods. At a sufficiently high shear rate, the fully nematic, 
homogeneous state is stable. The shear flow is then suddenly switched off, after which 
the system becomes either unstable or meta-stable. As a system we use suspensions of 
/^-viruses, which are mono-disperse and very long and thin somewhat flexible particles. 
The equilibrium phase behavior for these semi-flexible rods as far as the binodal points are 
concerned, is well understood on the basis of Onsager theory, extended to include charge 



and flexibility 



Polymer is added to the dispersion in order to widen the region 



of isotropic-nematic phase coexistence, which renders phase separation experiments feasible 
In a previous study we obtained the spinodal point that separates the unstable and 
meta-stable region relevant for the initial state in the experiments described in the present 
paper In the present paper, experiments are performed at concentrations such that the 
cessation of shear will render the system unstable. We interpret our data on the basis of the 
recent microscopic theory [6] mentioned above. Experimental evidence is given that phase 
separation is indeed driven by orientational ordering, which enslaves the concentration. 
We use confocal microscopy to confirm that demixing is indeed proceeding via spinodal 
decomposition, rather than nucleation and growth, through the observation of an initial 
interconnected structure of inhomogeneities. In addition we also use SALS experiments 
because these have a better time resolution and a better statistics. The present experiments 
are qualitative in the sense that only the specific wave vector dependence of the unstable 
eigen mode is discussed, without systematically varying the amount of added dextran which 
leads to attraction between the rods. 

This paper is organized as follows. First we give a brief overview of the microscopic theory 
on SD of rod dispersions (| , leading to predictions that will be tested experimentally. After 
the section on materials and methods, experimental results are presented for both techniques. 
In the discussion we analyze our results using the predictions of the microscopic theory. 
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II. THEORY 



The time-evolution of the probability density function (pdf) of the orientations and po- 
sitions of an assembly of N rods is described by the so-called Smoluchowski equation. From 
this microscopic equation of motion, an equation of motion for the number density p(r, u, t) 
of rods at r with orientation u at time t can be derived by integration 
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-p(r,u,t) = -DV-D(u)-{Vp(r,u,t)-/?p(r,u,t)F(r,u,t)} 

+ D r K- {^(r,u,f) -/3p(r,u,*)T(r,u,*)} . (1) 

Here, D and D r are the orient at ionaly averaged translational diffusion coefficient and rota- 
tional diffusion coefficient of a non-interacting rod, respectively. The orientational depen- 
dence of the translational diffusion coefficient of a non-interacting rod is described by the 
tensor D(u) = I + uu . Furthermore, TZj (•••)= X Vft (• • • ) is the rotational oper- 
ator with respect to Uj, where V^. is the gradient operator with respect to the cartesian 
coordinates of Uj. For very long and thin rods with hard-core interactions, the average force 
F and torque T on a rod with position r and orientation u due to interactions with other 
rods are given by, 

F(r,u,t) = -VV eS (r,u,t) , and , T(r,u,t) = -U V eS (r, u, t) , (2) 

where the "effective potential", 

V eS = | DL 2 /r 1 y du' |u x u | J dl j dl' p(r+lLlu+lLl'u',u',t) , (3) 

has been introduced earlier by Doi and Edwards 19]. 

In order to describe initial decomposition kinetics, the density p(r, u, t) is written as, 

p(r,u,t) = pP Q (vL,t) +5p(r,u,t) , (4) 

where 5p is the small deviation with respect to the initial probability density function 
p~P (u,t = 0), with p = N/V the average number density of rods. Note that Pq is gen- 
erally a function of time, which reflects the temporal evolution of alignment of the otherwise 
homogeneous system. Although we treat in this paper a quench from the nematic state, 
we will now proceed by assuming that the initial state is isotropic. Non-isotropic initial 
states require numerical analysis, since an appropriate (non-linear) equation of motion for 



Po should be solved simultaneously to the equation of motion for Sp . The general features 
of demixing are probably not very different for the different initial states. 

During the initial stage of demixing, Sp can be expanded up to second order in spherical 
harmonics as, 

Sp(r,u,t) = A (r,t) + A 2 (r,t) : uu . (5) 

The scalar A$ is proportional to the local number density of rods, while the tensor A 2 
describes the development of orientational order during demixing. As will turn out, the 
number density A is enslaved by the orientational contribution A 2 . 

Using these definitions in the equation of motion Eq.(jT|). the corresponding equations of 
motion for A and A 2 can be derived. These equations can be solved, leading to, 

MKt) = F —r T 1 r i r / L ^ kk : A 2 k, t , 

kk:A 2 (k,t) = kk:A 2 (k,t = 0) exp{-A (_) t} , (6) 

where (p is the volume fraction of rods, and L and D are their length and thickness, respec- 
tively. Furthermore, k is the unit vector along the wave vector k, and A^ - ) is the eigen value 
related to the unstable mode, which is equal to, 

A<"> = 6D r { 1 - \% v (1 - H| (kL) 2 ) } - <D(kL*) . (7) 

Note that for ^ (p > 4, this eigen value is negative for sufficiently small wave vectors, 
so that, according to eq.flHJ), inhomogeneities with the corresponding wave length 2ir/k will 
grow in time without any time delay. This concentration marks the location of the isotropic- 
to-nematic spinodal and is in accordance with Onsager's prediction [lj]. 

The proportionality of the density A to the orientational contribution A 2 in Eq.fjHJ) 
reflects the enslavement of density to orientational order during demixing. That is, the 
transition is driven by orientational diffusion rather than translational diffusion. 

Although Eq.(J2|) has been derived for an initial isotropic distribution, we think that the 
main conclusion, i.e. orientational fluctuations dominate the phase separation, is also valid 
for phase separation starting from the nematic state. 

In a scattering experiment the total scattered intensity is related to the quantities A and 
kk:A 2 as, 

A (k, t) (1 - i (kL) 2 ) - ^ (kL) 2 kk : A 2 (k, t) ( 1 - Jo (kL) 2 ) + O ({kLf) 1 ' .(8) 
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As for gas-liquid spinodal demixing suspensions of spheres, the scattered intensity during 
isotropic-nematic demixing suspensions of rods exhibits a ring-like pattern where a particular 
finite wavevevector grows most rapidly. The occurrence of a maximum in the scattered 
intensity at finite wavevectors during demixing has a fundamentally different origin for 
demixing rods as compared to spheres. For gas-liquid demixing of suspensions of spheres, 
the eigenvalue (here referred to as ') itself exhibits an extremum at a finite wavevector. 
For spheres A^ - -* is of the form D k 2 [1 — a k 2 ], where D and a are wavevector-independent, 



positive coefficients [20. |2l|. The prefactor k 2 signifies the fact that diffusion of spheres over 
long distances takes longer times, while a k 2 signifies the stabilization of large concentrations 
gradients. It is easily verified that an eigenvalue of this form exhibits an extremum at a 
finite wavevector. For gas-liquid demixing of suspensions of spheres this maximum in the 
growth rate results in the maximum in the scattering pattern. For the isotropic-nematic 
demixing of suspensions of rods, the eigenvalue is of the form D [1 — a k 2 } (see eq.(J7J), that 
is, the prefactor k 2 as compared to spheres is missing here. As a consequence A*- - ** remains 
finite at zero wavevectors. This difference in wavevector dependence of the eigenvalue A^ 
for gas-liquid and isotropic-nematic demixing is due to the fact that gas-liquid demixing 
is governed by translational diffusion while isotropic-nematic demixing is (predominantly) 
governed by rotational diffusion. The maximum in the scattering pattern is now due to the 
combination of the wavevector dependence of the eigenvalue and the wavevector dependent 
prefactors to the time-exponent. Note that according to eq.©, the prefactor of the time- 
exponent is indeed ~ k 2 , rendering the scattered intensity equal to zero at zero wavevector, 
which expresses conservation of the number of rods. The fc-dependence of X^ can be 
tested experimentally using the fact that, according to Eqs. (j6l8J) . the scattered intensity is 
~ exp{A^t}. Hence, 

| \n{I(k, t)} = 2 A" = 12 D r { 1 - I % <p (1 - H! (kL) 2 ) } . (9) 

The slope of a plot of ln{I(k,t)} as a function of t for a given wave vector is thus equal 
to A*"-* for that particular wave vector. Repeating this for various wave vectors allows to 
construct the wave vector dependence of X^\ 

In addition, the critical wave vector k c above which the system becomes stable, that is, 
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where ' becomes positive, is equal to, 




For shallow quenches, that is, for concentrations where the (p is close to 4, the critical 
wave vector is thus relatively small. That is, shallow quenches result in relatively large scale 
inhomogeneities while deeper quenches give rise to relatively small scale inhomogeneities. 

In this paper we prepare an initial nematic state, by shearing a suspension at large enough 
shear rate such that the induced nematic phase is stable against phase separation (see ref . |22j 
for a discussion of the bifurcation diagram for sheared systems), and then quench to zero 
shear rate. Since the orientation of the rods dominates the phase separation, it is expected 
that phase separation takes place anisotropically for an initial nematic state. In the following 
we will test this assumption and also the predictions made above for the isotropic-nematic 
SD, which we believe to hold true also for the nematic-isotropic SD. 



III. MATERIALS AND METHODS 

As model colloidal rods we use fd- virus particles which were grown as described in ref. Q| • 
The physical characteristics of the bacteriophage fd are: length L = 880 nm; diameter 
D = 6.6 nm; persistence length 2.2 \im. A homogeneous solution of 22.0 mg/mL fd- virus 
and 10.6 mg/mL of dextran (507 kd, Sigma-Aldrich, radius of gyration of 18 nm) in 20mM 
tris buffer at pH 8.15 with lOOmM NaCl is allowed to macroscopically phase separate. 
Without dextran the binodal concentrations are 21 and 23 mg/mL for the isotropic and 
nematic phase, respectively. Due to the added dextran, the binodal points shift to 17 and 
31 mg/mL, respectively. The lower spinodal point C spm for this sample was found to be 
equal to 24.7±1.1 mg/mL, as determined in a previous paper Q|. We prepared three 
dispersions by mixing a known volume of coexisting isotropic and nematic bulk phases of 
the quiescent dispersion. In this way, the osmotic pressure is independent of the varying 
ratio of dextran to fd- virus concentration. The mixing ratios of isotropic and nematic phases 
are chosen such that a quench from the nematic phase under flow will always render the 
aligned system unstable without flow, that is, the fd concentration is larger than the lower 
binodal concentration C bm = 17 mg/mL and smaller than the lower spinodal concentration 
(jspm = 24.7 ±1.1 mg/mL [18]. The fd concentrations are denoted here after as <pf, 
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where f = (C — C bm ) / (C spm — C bm ) relates to the fraction of the concentration between 
the lower binodal and spinodal. The concentrations used in the present study are : </>o.52 = 
19.3 mg/mL, 0o.55 — 19.9 mg/mL and 0o.84 — 23.6 mg/mL. For the SALS measurements 
the concentration of fd and dextran were 21.0 mg/mL and 12.1 mg/mL, respectively. Due 
to the fact that we used different concentrations of dextran for the different experiments, 
we cannot directly compare SALS data to the microscopy data. 

For the microscopy experiments we used a home-built counter rotating cone-plate shear 
cell, placed on top of a Leica TCS-SP2 inverted confocal microscope. This cell has a plane 
of zero velocity in which objects remain station ary with respect to the microscope while 
shearing. For details of the setup we refer to ref. 23j] . For the measurements described here 
we used confocal reflection mode at a wavelength of 488 nm. Quench experiments were done 
as follows. Samples were first sheared at a high rate of 10 s _1 for several minutes. The shear 
was then suddenly stopped, after which images were recorded at regular time intervals. For 
the SALS measurements we used a home made cylindrical optical shear cell. The rotating 
hollow inner cylinder has a radius of 21.5 mm, the gap width is 2.47 mm. The shear cell is 
placed in a cylindrical toluene bath with the second gap of the cell exactly in the middle of 
the bath. A 5mW He-Ne laser (Melles-Griott) with a wavelength of 632.8 nm was used as a 
light source. The laser beam is directed along the gradient direction through one single gap 
using a periscope system, which is inserted into a silicon oil filled inner cylinder at a fixed 
position. In this way the flow-vorticity plane is probed. Scattered intensities are projected 
on a white screen, with a beam-stop in the middle . The size of the beam-stop corresponds 
to a scattering angle of 1.4° and a wave vector of 2.4 x 10 5 m _1 . Images were taken in 
transmission with a peltier cooled 12 bit CCD camera, with 582x782 pixels (Princeton 
Instruments, microMAX). The maximum scattering angle was 5.1°, which corresponds to a 
wave vector of 1.0 x 10 6 m _1 . 



IV. EXPERIMENTAL RESULTS 

Confocal images of the morphology during spinodal decomposition in the early stage after 
a shear rate quench for different concentrations are given in the top row in FigH while in the 
bottom the corresponding Fourier transforms are plotted. The observed ring-like scattering 
patterns, typical for spinodal decomposition, are anisotropic with symmetry around the 
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[fcf]=23.6 mg/mL; 69 s [fcf]=19.8 mg/mL; 44 s [fc/]= 19.3 mg/mL; 31 s 




FIG. 1: (Color) The initial stages of SD for three different concentrations, as indicated in the 
figure. The top row shows the micrographs taken by reflection confocal scanning laser microscopy 
(field of view = 110 /im); the bottom row shows the fourier transform of the micrographs. 




FIG. 2: The scattering pattern for a sample with [fd]= 21.0 mg/mL and [dextran]= 12.1 mg/mL 
taken 23 s after cessation of flow. 
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nematic director which is along the flow direction. The same anisotropy is observed in the 
SALS measurements, where such Fourier space images are directly probed (see Fig[2J). 

Cross sections of the Fourier transforms and scattering patterns parallel and perpendicu- 
lar to the director are given in FigJ3]a and b, respectively. From these profiles we obtain the 
wave vector at which the Fourier transform exhibits its maximum, k max L, quantifying the 
wavelength of the fastest growing Fourier component of the inhomogeneous morphology. The 
values for k max L for directions parallel and perpendicular to the director as obtained from 
such cross sections of the Fourier transform of the confocal images and the SALS patterns 
are plotted as a function of time in FigEk and b, respectively. The ratio k maXt ±/k maXt \\ of 
k m ax perpendicular and parallel to the director is plotted in FigEt- This figure thus charac- 
terizes the anisotropy in the morphology. Since this ratio is larger than one, the typical size 
of inhomogeneities perpendicular to the director is smaller than the size of inhomogeneities 
in the direction parallel to the director. 

Comparing the profiles obtained from microscopy and from SALS, it is obvious that the 
SALS signal is less noisy. This is due to the fact that the volume that is being probed by 
microscopy is much smaller than the volume probed by SALS. On the other hand, due to the 
beam stop the profile starts at higher k-values than the profiles obtained from microscopy. 
The large difference in the scattered intensity parallel and perpendicular to the flow direc- 
tion is also observed under identical conditions for isotropic systems, and is caused by the 
experimental set-up ^J]- 

V. DISCUSSION 

The main result of the theoretical treatment in section |H] for SD of a dispersion of rods 
is expressed by Eq.fjHJ). This equation shows that changes in the concentration of rods are 
enslaved by changes in the orientation of the rods. Although the treatment was done for an 
initially isotropic state we believe that this result is valid independent of the initial condition. 
This relation between concentration and orientation has a few consequences, which can be 
tested experimentally: 

1) as for spheres there is no demixing for k — > 0; 

2) but the driving force of phase separation given by A*"-* goes to a constant value for 
k — > and not to zero as would have been the case when translational diffusion dominates 
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FIG. 3: Cross sections of the Fourier transform of the confocal images parallel and perpendicular 
to the director for a /^-concentration of ^0.52 (a), and the parallel and perpendicular cross sections 
of the scattered intensity as found from SALS (b). The wave vector is scaled by the length L of 
the fd- virus, the intensities are given in arbitrary units. 

phase separation kinetics (see the exponent in Eq.©). As a result k max shifts in time also 
in the initial stage of demixing; 

3) shallow quenches result in relatively large scale inhomogeneities while deeper quenches 
give rise to relatively small scale inhomogeneities, see Eq. fjlUj) . 

All micrographs in images FigQ show interconnected structures typical for spinodal de- 
composition. The resulting Fourier transforms in FigQas well as the SALS pattern in FigO 
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FIG. 4: The wavelength of the fastest growing Fourier component k max L as found from the cross 
sections parallel to the director (open symbols) and perpendicular to the director (filled symbols) 
of (a) the Fourier transforms of confocal images with a /^-concentration of 0o.84 (squares) and 0o.52 
(circles), and (b) the SALS pattern, (c) The anisotropy of the spinodal structure given by the ratio 
k m ax,±L/k max \\L for 0o.84 (squares) and the SALS measurement (triangles). 

are slightly elongated in the direction of the director. The observation of a ring structure 
confirms the prediction that the scattered intensity is zero for A; — 0, i.e. that there is no 
demixing for k — > (see also the cross sections in Fig|HJ). As pointed out in the theory 
section this is a general feature of SD, and it is a consequence of the conservation of the 
number of rods. The anisotropy in the morphology as well as in the growth rates, see Fig0t, 
show that the formation of inhomogeneities is affected by the initial orientation of the rods. 
The anisotropy in the phase separation is also seen in the case of nucleation-and-growth, 
where we observed that the nucleating factoids of isotropic phase are oriented along the 
director of the nematic background phase [18]. This anisotropy is due to residual alignment 
after the quench of the initially strongly sheared suspension. 

When plotting the wave vector where the intensity ring exhibits its maximum, i.e. k max , 
as a function of time, it is readily seen that this maximum shifts to smaller values also during 
the initial stage of demixing right after the quench (see Fig. 0]). That the initial stage of 
demixing is probed follows from FigEJ where the logarithm of the intensity is plotted versus 
time. In the initial stage this relation should be linear (see Eq.Q), which is indeed seen 
to be the case for all values of kL. The initial stage of the SD ends where this curve 
starts to deviate from linearity. For more shallow quenches closer to the spinodal, at higher 
concentrations, the initial stage extends up to 100 seconds. 

More importantly, for each value of kL we obtain the phase separation rate A* - -* from 
the slope of the curve (see Eq.(jOJ)). The resulting curves of A^ - ^ vs. kL are plotted in 
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FIG. 5: The logarithm of the Fourier component of confocal images for sample </>o.52 (a) and the 
scattered intensity from SALS experiments (b) as a function of time for various values of kL. The 
eigen value A^ - ) is obtained from the initial slope of this plot. 

FiglHlfor two different confocal microscopy samples (a) and the SALS sample (b). Clearly, 
A^ - ' approaches a finite value for k — > 0, as was predicted by theory, see FigJHb. This 
shows that the demixing kinetics is dominated by rotational diffusion. The absolute value 
for A^ - -* is about a factor ten higher for the SALS experiment than for the microscopy 
experiments, which is probably due to the lower dextran concentration that has been used 
in the microscopy experiment as compared to the SALS experiment. The analysis we used 

or our data could in principle also be applied to the measurements of van Bruggen et al. 

~y\ . where a similar behavior of the SALS patterns is observed. 

FigOJi also confirms the theoretical prediction in Eq. (|10jl that the initial inhomogeneities 
are larger for shallow quenches than for deep quenches. The length scale of the initially 
formed structures for the deep quench (circles, 0o.84) is H times the rod length in the direction 
parallel to the initial director and 7.1 times the rod length in the direction perpendicular 
to the initial director. For the shallow quench (squares, (po.52) the initial structure is barely 
anisotropic and has a typical size of 17 times the rod length. After about one minute the 
typical sizes for both concentrations start to overlap. 

Finally we would like to remind that the theory presented in section|n]is valid for repulsive 
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FIG. 6: The fcL-dependence of the eigen value A^ - -* as obtained from confocal microscopy for 
00.52 (bullets) and </>o.54 (squares) (a), as obtained from SALS (b), and as obtained from theory 
for the isotropic-nematic transition for different dimensionless concentrations jj(f)(c). Here L and 
D are the rod length and thickness, respectively, 4> is the volume fraction of the rods, and D r is its 
rotational diffusion at infinite dilution. 

rods. In our experimental system, however, depletion attractions between rods are induced 
by adding dextran. On adding more dextran, translational diffusion could play a more 
important role during phase separation. Rod-polymer mixtures has been treated on the basis 
of Ginzburg-Landau equations of motion with a thermodynamics input, amongst others, 
by Liu and Fredrickson Matsuyama et al.jnj, and by Fukuda 25]. In these papers 
it is shown that translational diffusion indeed becomes more important on increasing the 
polymer concentration. This results in a minimum in function of the wave vector at 

sufficiently high polymer concentration. We therefore believe that the pronounced minimum 
that is observed for the SALS measurement, see FigJHbj is related to attractions between the 
rods as induced by the added dextran. The two microscopy measurements depicted in FigEK 
are done at a lower overall dextran concentration. The sample with the somewhat lower 
fd concentration, i.e. higher polymer concentration, shows a less pronounced minimum as 
compared to the SALS sample, while for the sample with the lowest polymer concentration 
no minimum is present. These observations confirm the theoretical prediction on the polymer 
dependence of the phase separation kinetics. In future studies we will systematically vary 
the polymer concentration for a fixed concentration of fd. 
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VI. CONCLUSION 



We studied the nematic-isotropic SD of dispersions of fd-virus particles with added poly- 
mer after shear quenches into the two-phase region for varying concentrations. We exempli- 
fied the fundamental difference between spinodal decomposition of dispersions of rods and 
spheres using a recently developed theory. The main difference is that in the case of rods the 
phase separation is dominated by rotational diffusion. As a result the eigen value X^~\ which 
quantifies the rate of the phase separation, approaches a non-zero constant value for k — > 0, 
contrary to gas-liquid demixing of spheres where the corresponding eigen value becomes zero 
for k — > 0. This is due to the fact that for rods a local reorientation is sufficient to start the 
phase separation, whereas for spheres translational diffusion over finite distances is needed. 
We found experimentally the same fc-dependence of as predicted by theory 
Our experiments thus confirm that demixing is dominated by rotational diffusion and not 
by translational diffusion as suggested in earlier work j?l Q . In addition, we found a possible 
effect of translational diffusion through the minimum of the wave vector dependence of the 
unstable eigen mode, due to attractions between the rods as induced by the added dextran 
refs. ^, Q. This will be subject of further investigations. 
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